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Analysis of UBV Photometry in Selected Area 133 
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Abstract. We interpret the de-reddened UBV data for the field SA 133 to deduce the stellar density and metallicity distribu- 
tions function. The logarithmic local space density for giants, D*(0) = 6.40, and the agreement of the luminosity function 
for dwarfs and sub-giants with the one of Hipparcos confirms the empirical method used for their separation. The metallicity 
distribution for dwarfs gives anarrow peak at [Fe/H] = +0.13 dex, due to apparently bright limiting magnitude, V = 16.5, 
whereas late-type giants extending up to z ~ 4.5 kpc from the galactic plane have a multimodal distribution. The metallicity 
distribution for giants gives a steep gradient d[Fe/H]/dz — —0.75 dex kpc -1 for thin disk and thick disk whereas a smaller 
value for the halo, i.e. d[Fe/H]/dz — —0.45 dex kpc -1 . 
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1. Introduction 

Star count analyses of Galactic structure have provided a 
picture of the basic structural and stellar populations of the 
Galaxy. Examples and reviews of these analyses can be found 
in Bahcall (1986), Gilmore & Reid (1983), Majewski (1993), 
and very recently in Robin, Reyle, & Creze, (2000), Robin et 
al. (2003), Chen et al. (2001), and Siegelet al. (2002). Most of 
these programs have been based on photographic surveys; the 
Basel Halo Program (Becker 1965) has presented the largest 
systematic photometric survey of the Galaxy (Fenkart 1989a- 
d; Del Rio & Fenkart 1987; Fenkart & Karaali 1987;1990). 
The Basel Halo Program photometry is currently being recal- 
ibrated and reanalysed, using an improved calibration of the 
RGU photometric system which comprises homogeneous 
magnitudes and colours of about 20000 stars in a total of four- 
teen fields, distributed along the Galactic meridian through 
the Galactic centre and the Sun. While the ensemble of the 
full survey data are being analysed by comprehensive model- 
ing of the density, luminosity, and metallicity distributions of 
the different Galactic population components (Buser, Rong, 
& Karaali (BRK) 1998; 1999), the data in each individual 
field are being used as well for a detailed study of the lumi- 
nosity function (Ak, Karaali, & Buser 1998; Karata§, Karaali, 
& Buser 2001) and metallicity gradient (Rong, Buser, & 
Karaali 2001; Karaali et al. 2003; Karaali, Bilir, & Buser 
2004). In this paper we present the investigation of an individ- 
ual field in the U BV photometry, as some other works carried 
out as a part of Basel Halo Program, due to its Galactic posi- 
tion (I = 6°. 5). The abundance data have been derived from 
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broad band photometric estimates, primarily using the UV- 
excess (Sandage 1969). The technique of UV-excess to derive 
metal abundances has been applied for the photographic data 
in some works (e.g. Gilmore, Wyse, & Jones 1995; Karaali et 
al. 2004) though it broadens the metallicity distribution. 

Section 2 details the general characteristics of the photo- 
metric data set and the method. Section 3 discusses the den- 
sity function with two galactic models and comparison of the 
resulting luminosity function with that of Hipparcos (Jahreiss 
& Wielen (JW) 1997) and Gliese & Jahreiss (GJ 1992). In 
Section 4 we discuss the metallicity distribution. Section 5 
provides a conclusion. 

2. Data and Method 

2.1. De-reddening of Photometry 

The UBV data of 1729 stars in the field SA 133 (I = 6°. 5, 
b = +10°. 3) were taken from the Basel Photometric Cata- 
logue No. VIII (Becker et al. 1982) and the distances of 137 
stars, brighter than V = 15 mag, to the standard (U — B) a — 
(B — V) main sequence along the reddening line (Fig. 1) are 
used for reddening estimation. Among these stars 122 have 
only photographic al UBV data whereas the remaining 15 
stars have photoelectrical U BV data, taken from Mermilliod 
& Mermilliod (1994), additional to their photographical ones. 
The mean colour excess for all of these stars is E(B — V) — 
0.25 mag and their standard error for the mean, a = 0.01 
mag. This value is consistent with E(B — V) = 0.33 and 0.37 
mag derived for our field by Schlegel, Finkbeiner, & Davis 
(1998), and Burnstein & Heiles (1982), respectively. How- 
ever E(B — V) = 0.18 mag, given by Fenkart et al. (1986) is 
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Fig. 1. Two-colour diagram for stars used for estimation of 
the colour-excess, E(B — V) = 0.25 mag. 
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Fig. 2. De-reddened colour magnitude diagram for the field 
SA 133. 

less than our value. Thus, corrections applied to U — B and 
V due to reddening are E(U — B) = 0.72E(B -V) = 0.18 
and Av = 3.0E(B — V) = 0.75, respectively. The proce- 
dure applied here for de-reddening of the UBV data is the 
same as the one applied for all fields investigated in the Basel 
Halo Program (Hersperger 1973; Becker & Svopoulos; 1976 
Becker & Fang 1982; Becker & Hassan 1982). 

The catalogue error given by Fenkart et al. (1986) for U, 
B, and V magnitudes, for the magnitudes brighter than 17, is 
± 0.02 mag which corresponds to ± 0.03 mag in B — V and 
U — B. The colour magnitude diagram, V — (B — V) , in 
Fig. 2 indicates a limiting magnitude of V a — 16.5 mag for 
our field. 



with the direction of our field. The positions of some stars 
in our two-colour diagrams seem to be affected by the scat- 
tering which is unavoidable in the photographic photometry 
for faint stars. The most conspicuous appearance of these 
stars are those occupying the position of the most metal- 
rich stars though they are not so large in number to affect 
our results. Contrary to the metal-rich stars, extreme metal- 
poor stars ([Fe/H] < —3.0 dex) are sparse, especially in 
apparently bright V— intervals. Such stars exist in all two- 
colour diagrams in the U BV and RGU systems and some- 
times they are in an unexpected considerable amount of num- 
ber (Karatag et al. 2001). Hence, they are excluded from 
statistics without regarding their identity which can be bi- 
nary stars, extra-galactic objects or their position may be af- 
fected by relatively large photographic errors which is the 
case for faint magnitudes. The iso-metallicity lines, for the 
range — 3 < [Fe/H] < +1 dex, in Fig. 3 were adopted from 
Lejeune, Cuisinier, & Buser (1997). 

The separation of dwarfs and evolved stars (sub-giants or 
giants) was carried out to obtain a luminosity function consis- 
tent with the local luminosity function of nearby stars due to 
Gliese & Jahreiss (1992) and Jahreiss & Wielen (1997). The 
procedure of this separation is based on the fact that the local 
luminosity functions obtained for many fields indicated a sys- 
tematic excess of star counts relative to the luminosity func- 
tion of nearby stars for the fainter segment, i.e. M(G) > 6 
mag, and a deficit for the brighter segment, M(G) < 5 mag. 
The excess cited is due to the contamination of evolved stars 
and it is difficult to separate them from the dwarfs due to 
the lack of a spectral band sensitive to the surface gravity, 
log g, in the U BV and RGU photometric systems. However, 
it was shown in some recent works, if apparently bright stars 
with M(G) > 6 mag on the main sequence are removed to 
the category of evolved stars, both segments cited above con- 
verge to the local luminosity function of nearby stars (Karaali 
1992; Ak et al. 1998; Karatas et al. 2001; Karaali et al. 2004). 
Because this process causes decrease in the number of abso- 
lutely faint stars, M{G) > 6 mag, whereas it increases the 
number of bright stars. The apparent limiting magnitude of 
bright stars considered, lies within 15 < G < 16 mag. The 
best fit between the local luminosity function deduced in this 
work and the local luminosity function of nearby stars re- 
sulted, from iterations, when stars apparently brighter than 
V = 15.5 mag and absolutely fainter than M(V) = 5.5 
mag on the main sequence were assumed to be evolved stars. 
Thus, brighter absolute magnitudes were attributed to them 
and they were separated to giant - or sub-giant - category 
according to their magnitudes, i.e. M(V) < 2 mag and 
M(V) > 2 mag, respectively. 



2.3. Absolute magnitude and metallicity determination 



2.2. Separation of dwarfs and evolved stars 

The F-fractioned two-colour diagrams in Fig. 3 show that 
there is a considerable amount of metal-rich stars consistent 



The absolute magnitudes and metallicities are evaluated by 
means of two methods. For dwarfs with [Fe/H] > —1.75 
dex, we adopted the metallicity-luminosity calibration of 
Laird, Carney, & Latham (LCL 1988) which is given in terms 
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Fig. 3. The T^-fractioned two-colour diagrams: (a) Vo < 13.0, 
(b) (13.0 14.0], (c) (14.0 14.5], (d) (14.5 15.0], (e) (15.0 15.5], 
(f) (15.5 16.0], and (g) (16.0 16.5]. Symbols: (•): dwarfs, 
(o): sub-giants, (x): giants, and (+): not included into statis- 
tics due to its position in the two-colour diagram. The iso- 
metallicity lines, —3 < [Fe/H] < +1 dex, were adopted 
from the synthetic data of Lejuene et al. (1997). 



of the offset in absolute V magnitude from the Hyades main 
sequence, for which they obtained 

M{V) [Hyades) = 5.64(B — V) + 1.11 (1) 

The metallicity-dependent offset from this fiducial sequence 
that LCL derived ( AM{f ) is, for a star of given (B - V) 
colour and given UV-excess 5; 

AMy =[ 2,31 ~ 1 1 °^^~ ][— 0.6888<5+53.14i5 2 — 97.004i5 3 ] 

LCL state this calibration to be valid for 5 < 0.25, which cor- 
responds to [Fe/H] > —1.75 dex with Carney (1979) trans- 
formation of 8 into [Fe/H], as used by LCL and recently by 
Gilmoreetal. (1995). 

[Fe/H]= 0.11 - 2.90(5 - 18.68 8 2 (3) 

As already cited, the calibration of LCL is valid only for 
dwarfs with [Fe/H] > —1.75 dex. Hence another method 
was necessary for additional metal-poor dwarfs and evolved 
stars of any metallicity, for absolute magnitude and metallic - 
ity derivation. With a optimistic approach, the (U—B, B—V) 
two-colour diagram calibrated for iso-metallicity lines (— 3 < 
[Fe /H] < +1 dex) for dwarfs and evolved stars individu- 
ally, by means of the synthetic data of Lejeune et al. (1997), 
can provide metallicities with appropriate interpolation. For 
absolute magnitude derivation of giant stars, we calibrate 
M{V) — {B — V) a relation, using published data from of var- 
ious Galactic globular clusters to cover a large range of metal 
abundances. The globular clusters used are M92 ([Fe/H] = 
-2.24 dex), M5 (-1.40), 47 Tuc. (-0.71), and M67 (0.00), re- 
spectively, (Stetson & Harris 1988; Richer & Fahlman 1987; 
Hesser et al. 1987; Montgomery, Marschall, & Janes 1993). 

For the metallicity calibration of Carney (1979), an error 
of ± 0.02 mag in 8 yields an uncertainty of 0.3 dex. The ac- 
curacy of the calibration of the Basel Library Spectra (Leje- 
une et al. 1997) is in level of less than 0.05 mag, implying 
an uncertainty in derived [Fe/H] of up to 0.35 dex. The un- 
certainties in derived metallicities from two method are ap- 
proximately the same. We assumed that the same holds for 
late-type giants. The uncertainty ± 0.03 mag in (B — V) Q , 
cited above, corresponds to ± 0.2 mag in absolute magnitude 
M{V) derived by the colour-magnitude diagrams of four 
clusters. Combination of the uncertainities in V a and M(V) 
in the equation gives the distance r to the star, i.e. 

V - M(V) = Uog{r) - 5 (4) 

amounts to an uncertainty of 10%, in the distance estimate, 
which is less than the one (20%) cited by Gilmore et al. 
(1995). 

3. Density and luminosity functions 

Logarithmic space densities, D* = logD + 10 are evalu- 
ated for dwarfs and sub-giants for consecutive absolute mag- 
nitude intervals, i.e. (2 3], (3 4], (4 5], (5 6], (6 7], and 
(7 8] (Table 1), and for late-type giants (Table 2), where 
D(r) = iV/AVi,2, N: number of stars in the partial volume 
AV1.2 = (7r/180) 2 (D/3)(r| - rf), n, r 2 : limiting distances 
of AV1.2, □: 0.19 square-degree, apparent size of the field, 
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Table 1. The logarithmic space densities D* for dwarfs and sub-giants for all population types (distances in kpc, volumes in 
pc 3 . The symbols are defined in the text. Underlines indicate limiting distance of completeness. 
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Table 2. The logarithmic space densities for late-type giants 
(symbols as in Table 1). The standard deviation is s = ±0.39. 
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5.96 
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8.55 
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4.93 


10.00-15.85 


5.75 (7) 


13.55 


108 


4.27 


15.85-25.12 


2.29 (8) 


21.48 


58 


3.40 


25.12-39.81 


9.12(8) 


34.05 


40 


2.64 



r = [(rf + r 3 ,)/2] 1 / 3 centroid distance of the corresponding 
partial volume AVi.2- 

The density functions are compared with the galac- 
tic model of Gilmore & Wyse (GW, 1985), and Buser et 
al. (BRK, 1998; 1999) given in the form MogD(r) = 
logD(r, I, b) — logD(0, 1, b) versus r, where AlogD(r) is the 
logarithmic difference of the densities at distances r and at 
the Sun. Thus AlogD(r) = points the logarithmic space 
density at r = which is available for luminosity function 
determination. The comparison is carried out as explained 
in several studies of the Basel fields (Del Rio & Fenkart 
1987; Fenkart & Karaali 1987), i.e. by shifting the model 
curve perpendicular to the distance axis until the best fit to 
the histogram results at the centroid distances (Fig. 4a-g). 
Both models involve two disk components, thin and thick 
disk, and the de Vaucauleurs spheroid. The parameters for 
two models are given in Table 3. Although the two set of 
parameters (especially the scale heights) are different from 
each other, the model gradients (solid line for the model 
GW and dashed line for the model BRK) are close to each 
other in Fig. 4a-f, probably due to relatively short distances 
(r < 4 kpc). Actually, the last histogram section for giants 
(Fig. 4g), 25.12 < r < 39.81 kpc, show that the model 
gradient for BRK begins to diverge from the model gradient 
for GW at larger distances. There is adequate agreement be- 
tween GW and BRK models and the observed density func- 
tions for dwarfs and sub-giants within the limiting distances 
of completeness marked by horizontal thick lines in Table 1 
and the same agreement results, when one includes the lo- 



Table 3. Parameters for the models of GW and BRK. The 

symbols give: nc local space density, Hf. scale-height, he 
scale-length, c/a: axis ratio for the halo (scale heights and 
scale lengths in pc). 
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Parameter 


GW 


BRK 


thin disk 
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1 


1 




Hi 


300 


290 




hi 


4000 


4010 


thick disk 


n 2 /ni 


0.02 


0.06 




H 2 


1000 


910 




h 2 


4000 


4250 


halo 


n 3 /ni 


0.001 


0.0005 




c/a 


0.85 


0.84 



cal luminosity function. The luminosity function close to the 
Sun: <f*(M), i.e. the logarithmic space density for the stars 
with M ± m .5 at r = is the D* -value corresponding 
to the intersection of the model-curve with the ordinate axis 
of the histogram concerned. Fig. 5 shows the agreement be- 
tween the luminosity function resulting from comparison of 
our space density data with the GW model and the local lu- 
minosity functions of GJ and JW. 

However there is some difference between GW model 
and the observed density function for late type giants (Fig. 
4g), extending up to r « 40 kpc (z « 7 kpc distance to 
the Galactic plane). Actually there are 35 stars in excess in 
the last distance interval which cannot be adopted as dwarfs 
due to the agreement of the luminosity function with the ones 
cited above. It is interesting that the same stars do show an 
excess in the comparision between BRK model and the ob- 
served density function. Although the gradients for two mod- 
els do not fit in every distance interval, i.e. the gradient for 
the model of BRK is slight above the model gradient of GW 
at r w 16 kpc and it is below for the distance interval r > 26 
kpc, they present the same local space density. We do not 
have enough information for the exact identification of these 
stars in excess. They may be related to the fact that a num- 
ber of stars have been eliminated because they have been as- 
signed a too faint metallicity ([Fe/H] < —3 dex) or they 
may be extra-galactic objects. If the stars mentioned in the 
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Fig. 4. Logarithmic space density functions for dwarfs and 
sub-giants for six consecutive absolute magnitude intervals: 
(a) (2 3], (b) (3 4], (c) (4 5], (d) (5 6], (e) (6 7], (f) (7 8]; and 
for late-type giants (g). Thin solid and dashed lines represent 
model gradient of GW and BRK, respectively. The last his- 
togram section (25.12 < r < 39.81 kpc) with dashed line 
corresponds to the density for stars in excess. 




Fig. 5. Luminosity function for SA 133 resulting from com- 
parison of observed histograms with the model gradient GW 
(•), and confronted to the luminosity functions of GJ (0), and 
JW (H). 
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Fig. 6. Metallicity distribution for dwarfs and sub-giants (a), 
and late-type giants (b). Curves in (a) and (b) are the fitted 
gaussian distributions. 



first case were not eliminated the observed densities in the 
preceding distance intervals would be increased and there 
would not be appeared stars in excess in the last distance in- 
terval. In this case the local space density would be closer 
the the one of Gliese (1969). The second case seems also a 
solution for the problematic 35 stars, because the number of 
extra-galactic objects increases when one goes to faint mag- 
nitudes. Unfortunately the identification of the extra-galactic 
objects in this field could not be carried out by the University 
of Minnesota due to its low Galactic latitude (b = +10°. 3). 
If we omit these stars in excess, we obtain a logarithmic local 
space density D*(0) — 6.40 between those of Gliese (1969) 
and Fenkart (1989c), i.e. = 6.64 and D*(0) = 6.29, re- 
spectively. The errors for the density functions are given in 
Tables 1 and 2, and in Fig. 5 in standard deviations. 



4. Metallicity distribution 

A gaussian fit to the metallicity distribution for dwarfs 
and sub-giants (Fig. 6a) shows a high and narrow peak at 
[Fe/H] = +0.13 dex and a long metal-poor tail, though 
with less contribution to the total metallicity, expected for a 
low-latitude field to the galactic centre direction. Contrary to 
dwarfs and sub-giants, metallicity distribution for late-type 
giants which extend up to r = 25 kpc relative to the Sun 
or z = 4.5 kpc above the galactic plane thus forming a sub- 
sample of different population types is multi-modal (Fig. 6b). 
A gaussian fit to Fig. 6b gives two peaks at [Fe/H] = +0.12, 
—0.83 dex, and additionally a flat metal-poor tail extending 
down to [Fe/H] = -3 dex. 
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Each histogram in Fig 7a-c is fitted a gaussian curve with 
a mean < [Fe/H] >= -0.14, -0.78, and -1.10 dex. z- 
heights shown in the panels are < z >= 0.96, 1.72, and 2.20 
kpc, respectively. One can notice that a systematic shift from 
metal rich stars to the metal poor ones in Fig. 7a-c, where the 
mean metallicity as a function of z distance is displayed. 

Comparison of the mean metal abundances and < z > 
distances with any of these panels gives metallicity gradients 
as follows 

d[Fe/H]/dz = [(-0.78) - (-0.14)]/(1.72 - 0.96) = -0.84 dex kpc" 1 
d[Fe/H]/dz = [(-1.10) - (-0.14)]/(2.20 - 0.96) = -0.77 dex kpc" 1 
d[Fe/H]/dz = [(-1.10) - (-0.78)]/(2.20 - 1.72) = -0.67 dex kpc" 1 



The gaussian of Fig. 7d shows a peak at < [Fe/H] >= 
— 1.35 dex, with a < z >= 3 kpc, which reveal three metal- 
licity gradients close to each other but different than cited 
above, when compared with the data in each former panels. 
The mean of these metallicity gradients, i.e. d[Fe/H]/dz = 
—0.45 dex kpc -1 , give the indication that metallicity gradient 
is less steep in the outward of the Galaxy. 

The mean z-distance 1 kpc corresponding to a mean dis- 
tance 5.7 kpc in the Galactic plane lies in the dominant re- 
gion of thin disk, and the corresponding metallicity gradient, 
d[Fe/H]/dz = -0.84 dex kpc -1 , is close to the one given by 
Trefzger, Pel, & Gabi (1995) (d[Fe/H]/dz = -0.55 ± 0.10 
dex kpc -1 ) up to z = 0.9 kpc for a field investigated by means 
of Walraven VBLUW photometry. These authors claim a 
less steep gradient, -0.23 ± 0.04 dex kpc -1 , for z < 4 kpc, 
however. The metallicity gradient for thin disk, deduced from 
the investigation of seven fields by Rong et al. (2001), very 
recently, is d[Fe/H]/dz > — 0.6 dex kpc -1 , not contradict- 
ing with our results, within the limits of accuracy. Stars in 
the third and fourth panels (< z > = 2.2 and 3 kpc) in our 
work sample the galactic halo for which metallicity gradient 
is not controversial. Thus, the metallicity gradient for stars 
in the second panel (< z > = 1.72 kpc) to be discussed, 
for the region occupied by these stars is dominant by thick 
disk. Such a metallicity gradient contradicts some formation 
histories postulated for the formation of the (classical) thick 
disk. Until recent years this component of our Galaxy was 
assumed to have a mean metal abundance [Fe/H] ~ —0.60 
dex, a scale-height H < 1 kpc, and its space number density 
2 - 8 % of the thin disk in the solar neighbourhood. Addi- 
tionally, and more important, it was argued that the stars of 
thick disk were formed from a merger into the Galaxy (cf. 
Norris 1996 and references within), a formation mechanism 
unlikely to leave an abundance gradient. Some recent works 
suggested that the thick disk is more important component of 
the Galaxy, extending up to ~ 3 - 5 kpc from the galactic 
plane (Majewski 1993) with a metal-poor (Norris 1996) and 
a metal-rich tail (Carney 2000; Karaali et al. 2000; Karaali et 
al. 2004). Hence a revision of the formation scenario of the 
thick disk may be required. Furthermore, the works of Reid 
& Majewski (1993) and Chiba & Yoshii (1998) also suggest 
a metallicity gradient for the thick disk. 




-3 -2 -1 -0 1 

[Fe/H] 



Fig. 7. Metallicity distributions for late-type giants for four 
z-intervals which reveal a metallicity gradient into the direc- 
tion of the Galactic centre. Curves in each panel are the fitted 
gaussian distributions. 

5. Conclusion 

The agreement of the luminosity function resulting from the 
comparison of GW model with the logarithmic space density 
functions for stars with 2 < M(V) < 8 mag, and the log- 
arithmic space density function for late-type giants with the 
same model (and also with the model BRK) confirm the sep- 
aration of dwarfs and evolved stars in our field. 

There is a concentration at [Fe/H] = +0.13 dex and 
a long tail down to [Fe/H] ~ —3 dex for the metallicity 
distribution for dwarfs and sub-giants as expected for a low 
latitude field into the galactic centre direction with a short ap- 
parent limiting magnitude, i.e. V =16.5 mag. Whereas late- 
type giants which lie up to large distances relative to the 
Sun or up to large z-heights from the galactic plane have 
multi-modal metallicity distribution. However, the weighted 
mean metal-abundances of late-type giants in our sample is 
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< [Fe/H] >= —0.7 dex, rather close to the one claimed by 
Morrison & Harding (1993), i.e. —0.8 dex, who investigated 
only the K giants in two square degree field almost symmet- 
ric to ours relative to the galactic plane (/ = 350°, b = —10°) 
and to Harding's (1996) finding, < [Fe/H] >= -0.8 dex, 
for G and K giants. 

In this work, we observed metallicity gradient in each of 
the populations. This is different than the one usually ex- 
pected: the metallicity gradient may appear when there is a 
mixing of populations for which the relative proportion of 
the populations changes along the line of sight. In such case 
it results in a gradient which may not be present in each of 
the populations. However, in the present data there is a vari- 
able proportion of the populations when z is increasing, indi- 
cating a gradient in the populations. Additionally, it is rather 
steep in the regions dominated by thin and thick disks, i.e. 
d[Fe/H]/dz = -0.75 dex kpc -1 . As cited in Section 4, 
Reid & Majewski (1993) and Chiba & Yoshii (1998) also sug- 
gest a metallicity gradient for the thick disk. All these works 
would have implications on the discussion about the scenario 
of formation of the thick disk. However, we recognize that 
there are significant statistical uncertainties in our resuls, due 
to photographic data. It is also probable that the direction of 
our field (/ = 6°. 5) may play a specific role on our results. 
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